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Abstract. Average properties of general inhomogeneous cosmological 
models are discussed in the Newtonian framework. It is shown under 
which circumstances the average flow reduces to a member of the standard 
Friedmann-Lemaitre cosmologies. Possible choices of global boundary 
conditions of inhomogeneous cosmologies as well as consequences for the 
interpretation of cosmological parameters are put into perspective. 



1. The Standard Averaging Hypothesis 

The description of the Universe in terms of homogeneous solutions of Einstein's 
or Newton's equations for gravitationally interacting matter, respectively, is to 
be considered a historical remnant: Clearly, homogeneous cosmologies, in par- 
ticular the family of isotropic Friedmann-Lemaitre models, were motivated on 
observational grounds at the time of their innovation indicating a fairly homo- 
geneous distribution of matter around us, but were not justified on theoretical 
grounds. While the assumption of global isotropy has survived the observational 
experience until today's assessment of an extreme isotropy of the microwave 
background, the assumption of local isotropy around every point in space (which 
implies global homogeneity for analytical fields) is not adequate; the homoge- 
neous models have to be replaced by inhomogeneous cosmologies which account 
for the highly inhomogeneous matter distribution observed today. 

The homogeneous solutions are generally known to be linearly unstable to 
perturbations in the matter variables. Therefore, it cannot be expected that 
the standard models are reliable ones for the prediction of longtime behavior 
of universal dynamics. However, these models still enjoy application as "back- 
ground models"; empirically used parameters like the Hubble constant, the den- 
sity parameter and others are still thought to be determined by a homogeneous- 
isotropic solution, independent of spatial scale. Observations on scales which 
are assumed to be "statistically fair" are interpreted accordingly. This nai've de- 
scription of the global dynamical properties of the inhomogeneous Universe can 
only be justified, if Friedmann-Lemaitre cosmologies model the average dynam- 
ics correctly. The cosmological principle, in this context, may be considered an 
averaging hypothesis which asks for approval. 
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Examining present-day models of large-scale structure, either analytical 
approximations or numerical N-body simulations, we have to conclude that all 
models (that are actually used) are constructed such to obey the standard aver- 
aging hypothesis. A priori, however, there is no reason to believe this unless the 
average flow is calculated and the assumptions, which reduce the general average 
dynamics to a model which is a member of the Friedmann-Lemaitre cosmolo- 
gies, are specifled. Moreover, it has been pointed out several times, especially 
in the recent past (Shirokov & Fisher 1963, Ellis 1984, Futamase 1989) that 
the inherent nonlinearity of the equations (here always referred to the equations 
of General Relativity) manifestly contradicts the expectation that dynamical 
evolution and averaging commute, i.e., that the average of an evolved matter 
distribution agrees with the evolution of the average. This may be true for so- 
lutions of linearized equations or, in a statistical sense, for restricted ensembles 
such as for models based on the assumption of homogeneous-isotropic turbu- 
lence (Olson & Sachs 1973), but will not be in general the case as soon as we 
describe nonlinear stages of the evolution of matter variables, which is claimed 
to be covered by current structure formation models. 

Futamase (1989) (see also: Bildhauer 1990, Bildhauer & Futamase 1991a, b, 
Futamase 1995, Seljak & Hui 1995) proposes a general relativistic model to cal- 
culate the "backreaction effect" of local inhomogeneities on the global expansion 
law. In his model the effect is largest for strongly anisotropic gravitational col- 
lapse, inhomogeneities generally accelerate the expansion of the Universe, and 
they do it in a scale-dependent way. That this result has apparently not reached 
the community which builds and analyses models of large-scale structure may 
be explained by the ambiguity involved in averaging general relativistic models. 
Indeed, the averaging procedure is not uniquely determined, the metric enters 
as an additional dynamical variable (it may be averaged or, instead, deformed, 
as suggested by Carfora & Marzuoli 1988, Carfora et al. 1990). 

Since the main reason for the existence of a nonvanishing "backreaction" 
(i.e., the nonlinearity of the system of equations) is also present in the Newtonian 
framework, we can address the problem much easier there, and we will be able 
to build a precise notion of how to deal with it. However, even the Newtonian 
limit of averaged GR dynamics has not to coincide with the averaged Newtonian 
theory on scales which are thought to be fully covered by a Newtonian treatment. 

2. General Expansion Law in Newtonian Cosmology 

2.1. Time— evolution and averaging do not commute 

Let us consider a portion of the Universe D(t) with volume V(t). Henceforth, 
we concentrate on the expansion which we describe by the local expansion scalar 
= V - v. Introducing the scale-factor via the volume, ajj := V^^'^, we can write 
the spatial average of 9 on the domain D as (Buchert & Ehlers 1995): 

{0)d = U d^x0 = ^ = 3^ . (1) 
V Jd V (Id 

As written in (1) the spatial average can be calculated as a simple Euclidean vol- 
ume integral over the domain D, the main advantage of a Newtonian treatment. 
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The subscript D indicates that the averages (as well as the scale-factor) depend 
on morphological properties of the spatial domain such as content, shape and 
connectivity. 

We evaluate the evolution of the average in a tube of trajectories of fluid 
elements, i.e., we introduce a Lagrangian mapping : if = f(X,t) which sends 
fluid particles from their initial (Lagrangian) position X to their final (Eulerian) 
position X. We use the Jacobian of this mapping, d'^x = Jd'^X , / := det 
to transform spatial averages to volume integrals in Lagrangian space^: 
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y jD(to) 



Using the Lagrangian time-derivative, 
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commutation rule (Buchert & Ehlers 1995): 
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Equation (3) shows that the evolution of the average and the average over the 
evolved field do not commute, their difference being 
given by the nonlinear fiuctuation term on the r.h.s. . 




Figure 1. The mapping from Lagrangian to Eulerian space is in gen- 
eral not diffeomorphic; it may not even be a homeomorphism (Figure 
inspired by Carfora & Piotrkowska 1995). 



2.2. The generalized Friedmann equation 

Averaging Raychaudhuri's equation for the evolution of the expansion scalar, the 
scale factor a£) is found to obey the general expansion law (Buchert & Ehlers 
1995): 

3^ + 4vrG^-A = -Q , (4a) 
an 

where the source term Q depends on the fiuctuation term in (3) and the mag- 
nitudes of rotation (uj) and shear (a) of the fiow, 

Q:=l{{O)l-{0')D)+2{a'-i.')D . (46) 



disclaimer is added in Figure 1. 
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Md denotes the total mass contained in D. 

Eq. (4) may be rewritten as a standard Friedmann equation for the actual source 
term {geff)D- 

4vrG'(^e//>D ■.= 4TrG{g)D + Q , (5) 

where {q)d = Mo/a'j) is the pure average matter density. Eq. (5) shows that, for 
irrotational flows (which we may consider a good assumption until the epoch of 
structure formation), the additional "dynamical mass" is a positive term which 
adds to the matter density, if is larger than the fluctuation (^)|) — = 
{{6 — {0)Dy)D > 0. This suggests to add the source term Q to the list of dark 
matter candidates: strongly sheared inhomogeneities could "fake" an additional 
density which, e.g., leads to an overestimate of the density parameter. 

Integrating eq. (4) with respect to time yields the generalized form of 
Friedmann's differential equation: 



^5 5 TT = T- / Q -TT^n • (6) 

a% 3a% 3 3al, Jt, ^ dt' ^ ^ ^ 



2.3. Averaging globally homogeneous— isotropic universes 

We now assume that a global Hubble flow exists on some large scale A and that 
the expansion-factor on that scale obeys Friedmann's differential equation (Eq. 
(6) for Q = 0); on the scale A we write a£) =: a. Splitting the velocity gradient 
Vij into its Hubble part and a peculiar-velocity gradient, Vij = H(t)Sij + Uij, 
where H(t) = ^, we obtain: 

= 3H + V -u . (7) 

After averaging and using Gaufi's theorem, the last equation leads to a relation 
between the Hubble function H(t), the "effective Hubble function" H£)(t) := 
and the peculiar-velocity field u(x,t): 

Hit) = Hoit) - \a-^^ j d^S-u . (8) 

3 JdD(t) 

Ho may be interpreted as that Hubble function which is inferred from the 
(possibly anisotropic and rotational) dynamics of the spatial domain D. (This 
interpretation is possible if statistical averages of many such spatial domains are 
considered, but, at present, we only measure one member of such an ensemble.) 

Accordingly, the source term Q can be split into its Hubble part and devi- 
ations thereof and transformed into surface integrals over the boundary dD{t) 
(Buchert k Ehlers 1995): 

3— + AttG^ -A=--(a7?[ d'S ■ u] + a7? [ dS- (uV -u-u-Vu) . 

an al 3 V ^ isD ) ^ JdD ' 

(9) 
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2.4. Fair samples and cosmological parameters 

In view of Eqs. (6), (8) and (9) the measured average characteristics like Hub- 
ble's constant, the density parameter, and others will, in general, show a scale 
dependence. The Hubble constant may be large on small scales where the effect 
of Q dominates, but small on larger scales; the density parameter shows the 
same tendency in conformity with the quantitative calculations by Bildhauer 
& Futamase (1991a, b). Consequently, the age of the Universe is likely to be 
underestimated on scales which we currently access observationally. Thus, look- 
ing at the presently discussed observational values of these parameters and the 
problems of large-scale structure models, this effect works in the right direc- 
tion: Models for large-scale structure favor a small Hubble constant, possibly 
as small as 30 km/sMpc as championed by Bartlett et al. (1995), observations 
of the local Universe point instead to a value of about 60 — 70 km/sMpc (Ellis 
et al. 1995), which has to be identified with the "effective Hubble function" as 
inferred from one spatial domain (our environment). 

Let us have a closer look at the "backreaction terms" in Eq. (9). Both 
terms are entirely neglected in the traditional approach, because one assumes 
that Jj-, d'^x V ■ u = on any postulated "fair sample" of the Universe (even if 
the sample is of the order of 100 Mpc shallow); also the other fiuctuation term is 
neglected by construction due to the assumption that the matter variables aver- 
age to zero on their assumed Friedmann backgrounds. Concerning the first term 
it is interesting to note that reconstruction models which try to recover the den- 
sity distribution from measurements of peculiar-velocities are determined only 
up to a constant offset factor (the mean density contrast {S)d on the observed 
scale, which is related to the square root of the first term in (9) using linear 
perturbation theory, S oc —V • u; see, e.g., Dekel 1994). 

The qualitative considerations above, however, indicate that estimates on 
intermediate and small scales will give the worst results for the mean character- 
istics of the Universe. The source term Q will be most effective in A— dominated 
universes where fiuctuations grow fast. To establish this effect quantitatively we 
could run large (Gpc-) simulations and look at sampling variations in the quan- 
tity [i {{0)jj - {0'^)d) + (o"^ - LL>'^)D]/'^T^G{g)D- Previous estimates of fiuctua- 
tions in the cosmological parameters were based on spherical symmetric models 
((7 = 0) which can only account for the fiuctuation term for 9 (e.g., Suto et 
al. 1995). Since structures form in a highly anisotropic way, the shear term 
will play an important, if not dominant role. It should be emphasized that not 
only amplitude effects matter which are mirrored in the scale-dependence of 
the r.m.s. density contrast or the power spectrum, but phase correlations can 
produce large (albeit low-amplitude) structures (Beacom et al. 1991, Buchert 
& Martinez 1993). As long as we do not cover scales which are considerably 
larger than the extent of these structures and as long as the phases are not 
uncorrelated, we can neither expect the "backreaction effect" to vanish, nor can 
we speak about a "fair sample" of the Universe. Generally speaking, since there 
is little parameter space left for the standard model (even in view of a conser- 
vative interpretation of observational and theoretical constraints, Bagla et al. 
1995, Ellis et al. 1995), there is a problem which could be due to the ignorance 
of "backreaction effects". 
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3. Global Interpretations and Alternative Views 



3.1. Toroidal cosmologies 

We may assume that, on the largest scale A, the cosmology is topologically 
closed, i.e. space sections are compact without boundary. The simplest example 
of such a space form is provided by the familiar hypertorus. As a consequence of 
this assumption the surface integrals in eq. (8) and (9) vanish and H£)(t) = H{t) 
on the scale A. (Note that the existence of a global Hubble flow represents a 
restriction of generality in the case of spatially compact models, because there 
might exist a global vorticity and/or shear flow on the scale A; see Ehlers & 
Buchert 1995.) Since toroidal cosmologies are realized if the perturbation fields 
are periodic on the box size of a simulation, we conclude that, technically, cur- 
rent structure formation models obey the standard averaging hypothesis for 
globally homogeneous-isotropic universes, i.e., they average out to the standard 
background models, which has been assumed so far without proof. However, 
the scale A is usually set ad hoc; the present analysis shows that it will be only 
meaningful to "close the topology" on a scale where the "backreaction terms" 
discussed have negligible infiuence, unless we deal with "small universes" (Ellis 
1971) having literally a closed topology. 

We have to be careful by extending this picture to General Relativity, e.g., 
for i7 < 1 universes the space sections have constant negative curvature and 
cannot be compactified to a hypertorus. (For details on possible space forms 
see, e.g., Ellis 1971, Lachieze-Rey & Luminet 1995 and ref. therein.) 

3.2. Infinite and topologically open cosmologies 

Although compact and unbounded cosmologies are attractive, if not favourable, 
it is possible that there doesn't exist any scale on which the "backreaction" 
vanishes, and there may be no Hubble fiow which fixes a standard frame of 
reference. However, we believe that the source term Q has smaller values on 
larger scales which may entitle us to "approximate" the Universe by standard 
models on the largest scales. However, the bare existence of Q indicates the 
presence of an intimate relationship between the evolution of inhomogeneities 
and the morphology of patches of smoothed-out matter distributions on every 
scale. To ignore this term entails an attitude towards "closing the eyes" in front 
of a possibly exciting development of inhomogeneous models of the Universe. 
A reinvestigation of globally hierarchical cosmologies is only possible with the 
source term Q. 

3.3. Scaling, coarse— graining and renormalization 

Armed with tools for a scale-dependent average dynamics, we might accom- 
plish a fully scale-dependent dynamical description of the Universe by advanc- 
ing renormalization group techniques. Having emerged from different branches 
of physics, these techniques might initialize the future of inhomogeneous cosmol- 
ogy, calling for the establishment of the "coarse graining" idea, which could shed 
light on general scaling properties of both gravitational dynamics and statistical 
characteristics of the matter distribution. The urgency of such an approach has 
already invoked various efforts (e.g., Carfora & Piotrkowska 1995, Piotrkowska 
1995, Perez-Mercader et al. 1995). 
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